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Abstract 

The present results on the Vr oscillation search by the CHORUS experiment at CERN are summarized. A 

fraction of the neutrino interactions collected in 1994-1997 has been analysed, searching for charged current 
interactions followed by the r lepton decay into a negative hadron, electron, or into a muon. A sample of 126,229 
events with an identified muon in the final state and 19,436 events without an identified muon in the final state 
have been located in the emulsion target. Within the applied cuts, no candidate has been found. This result 
leads to a 90% C.L. limit P{vf^ v-r) < 4.0 ■ lO""* on the mixing probability. Ave — >■ Vt exclusion plot is 
also presented corresponding to 90% C.L. limit on the mixing probability of P(!/e ^ i't) < 3.0- 10~^. 
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The CHORUS experiment has recently reported ^ 
I |] a Hmit on Vr oscillation obtained by the 

analysis of a subsample of neutrino interactions, taken 
in 1994-1996, both with and without an identified /i^ 
in the final state. This paper contains an update of the 
statistics of the subsample of events with an identified 
^~ in the final state. 

The experimental setup and the characteristics of 
the CERN wide band neutrino beam are summarized 
in ijlj and described in more detail in [^. 

1 The apparatus 

The " hybrid" CHORUS apparatus combines a 770 
kg nuclear emulsion target with various electronic de- 
tectors : a scintillating fibre tracker system, trigger ho- 
doscopes, an air-core magnet, a lead/scintillator calorime- 
ter and a muon spectrometer. 

Details about the experimental setup and the perfor- 
mances of the sub-detectors can be found in [^. 

2 Data collection and event selection 

In the 1994-1997 period, CHORUS has collected 
2,271,000 triggers corresponding to 5.06 10^^ protons 
on target. Of these, 458,601 have a muon identified in 
the final state (the so-called 1/i events) and 116,049 do 
not (the so-called 0/x events) and a vertex position com- 
patible with one of the four emulsion target stacks. 

All tracks, associated to the interaction vertex, with 
an angle smaller than 0.4 rad from the beam axis and 
bigger than 0.05 rad from the direction of a secondary 
beam nearby are selected for further analysis. The sec- 
ond requirement avoids possible confusion arising from 
the large emulsion background of muons originating from 
this beam. The selected tracks are searched for in the 
emulsion if their charge is negative and their momen- 
tum is in the range 1 < Pft, < 20 GeV/c and < < 
30 GeV/c for hadrons and muons, respectively. 

It should be noted that the energy deposition in the 
calorimeter has not been used to select electrons or uniden- 
tified muons. However the small contribution of the r~ 
leptonic decay modes to the Ofi data sample is taken into 
account in the evaluation of the sensitivity. 

3 Scanning procedure 
3.1 Vertex location 

The various steps leading to the plate containing 
the vertex by means of fully automatic microscopes are 
identical to those described in [|l| §]. They are applied 
to the muon for the 1/i events and to all the negative 
tracks for the 0/i events. Each track is searched for in the 
interface emulsion sheets positioned between the scin- 
tillating fibre tracker detector and the target emulsion 
stack. The found tracks are then followed upstream in 
the target emulsion, using track segments reconstructed 
in the most upstream 100 fim of each plate, until the 
track disappears. This plate is referred to as the vertex 
plate, since it should contain the primary neutrino ver- 
tex or the secondary (decay) vertex from which the track 
originates. The three most downstream plates of each 
stack are used to validate the matching with the inter- 
face emulsion sheets and are not considered as possible 



vertex plates. The scanning results are summarized in 
Table |. The mean efficiency of this scan-back proce- 
dure is found to be 30% for 0/x events and ^ 41% for 
1/i events. A detailed simulation of the scanning crite- 
ria shows that the difference mainly reflects the poorer 
quality of the hadron track predictions at the interface 
emulsion sheets, due to the difficulty of reconstruction 
inside a dense hadronic shower and the larger multiple 
scattering owing to the lower average momentum. 

3.2 Decay search 

Once the vertex plate is defined, automatic micro- 
scope measurements are performed to select the events 
potentially containing a decay topology (kink). Differ- 
ent algorithms have been applied, as a result of the pro- 
gress in the scanning procedures and of the improving 
performance in speed of the scanning devices. They are 
described in [|l], H and briefly recalled here. 

In the first procedure, the event is selected either 
when the scan-back track has a significant impact pa- 
rameter with respect to the other predicted tracks or when 
the change in the scan-back track direction between the 
vertex plate and the exit from the emulsion corresponds 
to an apparent transverse momentum, pr, larger than 
250 MeV/c. For the selected events and for those with 
only one predicted track, digital images of the vertex 
plate are recorded and are analysed off-line for the pres- 
ence of a kink. 

The second procedure is restricted to the search of 
long decay paths. In this case the vertex plate is assumed 
to contain the decay vertex of a charged parent particle 
produced in a more upstream plate. With this procedure 
only kink angles larger than 0.025 rad are detected. 

For the events selected by either one of these pro- 
cedures, a computer assisted eye-scan is performed to 
assess the presence of a secondary vertex and measure 
accurately its topology. A decay candidate must sat- 
isfy the following criteria: 

1. the secondary vertex appears as a kink without 
black prongs, nuclear recoils, blobs or Auger elec- 
trons; 

2. the transverse momentum of the decay muon (ha- 
dron) with respect to the parent direction is larger 
than 250 MeV/c (to ehminate decays of strange 
particles); 

3. the kink, in the 0/i channel, occurs within three 
plates downstream from the neutrino interaction 
vertex plate. Because of the lower background, 
the kink search in the muonic decay channel was 
extended to five plates, with a gain in efficiency 
of about 8%. 

No decay candidate has been found satisfying 
the selection criteria. 

4 Experimental check of the kink finding 
efficiency 

The kink finding efficiency, ekink, has been eva- 
luated by Monte Carlo simulation and experimentally 
checked looking at hadron interactions and dimuon events 



Table 1: Current status of the CHORUS analysis 







1994 


1995 


1996 


1997 


Protons on target 




0.81 • IQi'J 


1.20 • lO^'J 


1.38- IQi" 


1.67-101" 


Emulsion triggers 




422,000 


547,000 


617,000 


719,000 


1/i to be scanned 




66,911 


110,916 


129,669 


151,105 


0/i to be scanned 




17,731 


27,841 


32,548 


37,929 


1 /i scanned so far 




46,169 


52,130 


98,548 


108,796 


0/i scanned so far 




10,639 


13,364 


23,760 


19,344 


1/i vertex located (in 


the 33 most upstream plates) 


19,581 


21,809 


38,919 


45,920 


0/i vertex located (in 


the 33 most upstream plates) 


3,491 


4,023 


6,758^) 


5,1641) 



On a sample of about 55 m of hadron tracks scanned 
in emulsion 2 1 neutrino interaction events with a hadron 
interaction have been detected within the applied cuts 
during the decay search procedure. This result is in good 
agreement with the expected value of (24 ± 2) from 
Monte Carlo simulation. 

Part of the dimuon sample (two muons in the fi- 
nal state) collected in 1995 and 1996 has been analysed 
searching for the reaction I'^N pi~D^X with the 
subsequent muonic decay of the Assuming a charm 
yield of about ^ 5%, in the sample scanned we expect 
(2 2. 8 ±3. 9) dimuon events and we found 25 events. This 
result shows that Monte Carlo simulation and real data 
are in good agreement. Although the number of events 
is too low to draw quantitative conclusions, we can take 
these results as qualitative checks of the simulation of 
the automatic scanning procedure. 

5 Sensitivity and backgrounds 

In this section we discuss the expected background 
from known sources and, in absence of r~ candidates, 
we derive limits on the t/^ Vt and — > os- 
cillation parameters. Both signal efficiencies and back- 
ground estimation have been evaluated from large sam- 
ples of events, generated according to the relevant pro- 
cesses and passed through a GEANT [|| based simu- 
lation of the detector response. The output is then pro- 
cessed through the same reconstruction chain as used for 
data. The simulated tracks in emulsion are used to esti- 
mate the efficiencies of each scanning step. Apart from 
the kink detection efficiency, only ratios of acceptances 
enter the calculation of the experimental result and most 
of the systematic uncertainties of the simulation cancel 
out. 

5.1 Background estimates 

The main source of potential background in the muonic 
T channel is the charm production. We expect less than 
0.2 events in the current sample from the anti-neutrino 
components of the beam: 

For 1996 and 1997 the decay search on 0/i is not yet finished and 
the relative statistics are not included in the cuiTent limit 



i^M(^e)iV ^ n+{e+)D~X 
followed by 

D- n-X° 

in which the /i+(e+) escapes the detection or is not 
identified. 

The sources of background for the hadronic r decay 
channel are: 

- the production of negative charmed particles from 
the anti-neutrino components of the beam. These 
events constitute a background if the primary /i+ 
or e+ remains unidentified. Taking into account 
the appropriate cross-sections and the branching 
ratios, in the present sample we expect ^ 0.02 
events from these sources; 

- the production of positive charmed mesons in char- 
ged current interactions, if the primary lepton is 
not identified and the charge of the charmed par- 
ticle daughter is incorrectly measured. We expect 
~ 0.03 events from this source in the present sam- 
ple; 

- the associated charm production both in charged 
(when the primary muon is lost) and neutral cur- 
rent interactions, when one of the charmed parti- 
cles is not detected. The cross-section for charm- 
anticharm pair production in neutral current inter- 
actions relative to the total charged current cross- 
section has been measured [^] to be 0.131q'5i%. 
The production rate of associated charm in char- 
ged current interactions is unknown, but an upper 
Hmit is available (< 0.12% [^). In the present 
sample, the estimated background from this pro- 
cess represents < 0.01 events; 

- the main background to the hadronic t~ decays is 
due to so-called hadronic "white kinks", defined 
as 1 -prong nuclear interactions with no heavily 
ionizing tracks {black and grey tracks in emulsion 
terminology) and no evidence for nuclear break 
up (evaporation tracks, recoils, blobs or Auger elec- 
trons). Published data allowing to determine the 
white kink interaction cross-section are scarce wi 
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1^]. In a dedicated experiment with 4 GeV pions 
at KEK a very steep fall-off in pr was ob- 
served and only 1 out of 58 observed kinks had a 
PT larger than 300 MeV/c. Since the experimen- 
tal information of the pr dependence is statisti- 
cally poor at large values, a Monte Carlo simula- 
tion, based on a modified version of FLUKA 
[l0| ], has been performed. The results of this sim- 
ulation are in good agreement with the pt de- 
pendence of the KEK measurement. An effective 
white kink mean free path in emulsion, X^jk ^ 22 
m, is obtained for a px cut at 250 MeV/c, using 
the pion energy spectrum as observed in the 0/i 
sample. The above result is compatible with the 
4 observed events with pr > 250 MeV/c, at all 
distances from the primary vertex, along a total 
path of ~ 92 m of scan-back tracks. This cor- 
responds to a background estimate of 0.5 events 
within three plates downstream from the primary 
vertex for the 0/i sample of 1994 and 1995. The 
analysis on the larger statistics now available is in 
progress. From the study of the events with long 
apparent decay length we will be able to deter- 
mine, in a way largely free of systematics, the 
number of expected events with a kink within three 
plates from the primary vertex. In addition, the 
different kinematical properties of this background 
with respect to the r signal can be used, if neces- 
sary, to discard these events without spoiUng the 
sensitivity of the experiment. 
The prompt 1/7- contamination of the beam [ jlTI ] is 
a background common to both the hadronic and muonic 
decay channels. For the present sample the expected back- 
ground is less than 0. 1 events. 

5.2 Sensitivity to the i^^ Vr oscillation 

In the approximation of a two-flavour mixing scheme, 
the probability of appearance in an initially pure 
beam can be expressed as 



1) T ^ /i, 2) T — ^ ft,, 3) T — > e and 4) t — ^ /l (the fi 
is not identified) channels. 

The expected number, Nri (i — 1, 2, 3, 4), of ob- 
served decays into a channel of branching ratio BRi 
is then given by 



Nr,. = BR, • / $1.^ • P^,T ■ Or ■ Ari ■ Er* ' dE (1) 



sm 



2 / 1.27 ■ Aml^ieV^) ■ L{km) 
[ E{GeV) 



I'rVfj.n ) = (17.35 ± 




with 

- Si?(i„,4) = BR{t 

0.10)% [|12||); 

- BR2 = BR{t Vrh-nh^) = (49.78±0.17)% 

- Bi?3 = BR{t VrVee-) = (17.83±0.08)% ^ 

- $1/^ the incident flux spectrum; 

- (Ti- the charged current Vt- interaction cross-section; 

- Ari the acceptance and reconstruction efficiency 
for the considered channel (up to the vertex plate 
location); 

- tri the corresponding efficiency of the decay search 
procedure. 

With proper averaging (denoted by ()), Nri can also 
be written as a function of n, : 



Nr, = BR, ■ n, ■ {Pf,r) 



(gr) (Ari) 

(o,) ■ {A,) 



(2) 



where 



71 1 = (the number of located charged current 
interactions corresponding to the considered 
event sample) and 712 = 713 = 714 = (-/V^)o^ 
(the product of A^^ and the relative fraction of the 
0/i sample for which the analysis has been com- 
pleted); 

'Ofj.{T)) = J -di?. It takes into account 



fj.{Ti) 



dE 



dL 



where 

- E' is the incident neutrino energy; 

- L is the neutrino flight length to the detector; 

- Of^ir is the effective — Vr mixing angle; 

- Arn?^^ is the difference of the squared masses of 
the two assumed mass eigenstates; 

- '^{E, L) is the fraction of with energy E orig- 
inating at a distance between L and L + dL from 
the emulsion target. 



quasi-elastic interactions, resonance production and 
deep inelastic interactions iS{^^^)syst 

iSi^hyst - 7%); 
- (eri) is the average efficiency of the decay search 
procedure for the accepted events (S{{eri))syst ^ 
10%). 

To allow an easy combination of the results from 
the 1/i and O/i event samples, it is useful to define the 
"equivalent number of muonic events" of the 0/i sample 

as 



4 

E 

1=2 



(^r.) (£r.) BR, 

{Ar^,) ■ (e,^) ■ BR^ 



(3) 



The 90% C.L. upper limit on the oscillation proba- 
bility then simplifies to 



The T channels considered in the i/^ 
lation search we describe in this paper are: 



oscil- 



P, 



AIT 



< 



2.38 ■ 



' rA 



BR^ ■ (e,^) • [N^ + N^"] 



(4) 
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where r„ = (cr^)/(o-^) and = 
In the above formula, the numerical factor 2.38 takes 
into account the total systematic error (17%) following 



the prescription given in |13|. The systematic error is 



mainly due to the uncertainty of the Monte Carlo simu- 
lation of the scanning procedures. 

The estimated values of the quantities appearing in 
this expression are given in Table|. No statistical errors 
are quoted since they are much smaller than the system- 
atic uncertainty. 

Using the present sample the following 90% C.L. 
limit is obtained 



Pur < 4.0 ■ 10^ 



(5) 



In a two flavour mixing scheme, the 90% C.L. ex- 
cluded region in the {siTi?20^r, A^i^r) parameter space 
is shown in Figure [l|. Maximum mixing between and 
t/^ is excluded at 90% C.L. for A TO^^ > 0.7 eV^; the 
large Am^ are excluded at 90% C.L. for sm^26'pr > 
8 • 10-4. 

5.3 Sensitivity to the Vt oscillation 

The SPS neutrino beam contains a v^. component 
which amounts to 0.9% of the integrated flux. The 
negative result of the search for Vr interactions can there- 
fore be used to set limits on the Vf. Vt oscillation. 
The evaluation of the limit has been performed with the 
technique described in section 5.2, with the replacement 
of the energy dependent i^^ flux by the flux. The dif- 
ference between the j/^ and Vf, energy spectra - the aver- 
age energy of component is 20 GeV higher - leads to 
differences in the acceptance for interactions. In the 



case of the beam, the increase of the cross-section 
with energy improves the sensitivity to Vr interactions, 
while the kinematic cuts and reconstruction inefficien- 
cies affecting high energy events contribute to lower the 
acceptance. 

Using the present sample the foUowing 90% C.L. 
limit is obtained 



P,.r < 3.0 • 10" 



(6) 



The 90% C.L. excluded region in the (siTi^29eT, 
Aml^) parameter space is shown in Figure |[ Maximum 
mixing between I'e and i^r is excluded at 90% C.L. for 
A ml^ > 8 eV^; the large Am"^ ai-e excluded at 90% 
C.L. for sin^2eer > 6 • lO'^. 

6 Conclusions 

The emulsion scanning methods, previously described 
in have been applied to a fraction of the 1994- 

1997 data. No decay candidate has been found, lead- 
ing to a more stringent 90% C.L. upper limit on the 
f/^ Vt oscillation probability (P^r < 4.0 • lO"*). 
Based on the same statistics, the 90% C.L. upper limit 
on the Vf. v-r oscillation probability is Pf,T < 3.0 ■ 

10-2. 

The dominant background to the r search is due to 
"white kink" secondary interactions in the O/i channel. 
A direct measurement of this background process is in 
progress and kinematic cuts are planned for its reduc- 
tion. Furthermore, a second phase of the analysis (with 
better efficiencies, larger statistics and faster automatic 
emulsion scanning) has started with the aim of reaching 
the design sensitivity (P^,- < 1.0 • 10~*) 
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Table 2: Quantities used in the estimation of the sensitivity 





1994 


1995 


1996 


1997 


N 


19,581 


21,809 


38,919 


45,920 


T 

' a 


1.89 


1.89 


1.89 


1.89 


' A. 


0.93 


0.93 


0.93 


0.93 


\ 'H'l 


0.39 


0.39 


0.39 


0.39 


{Arh) 

\ • "-I 


0.17 


0.17 






{Are) 


0.093 


0.093 






{Ar-,) 


0.026 


0.026 








0.53 


0.35 


0.37 


0.37 




0.24 


0.25 






(ere) 


0.12 


0.13 








0.22 


0.23 








11,987 


14,769 
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